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ABSTRACT 

We present a survey of atomic carbon (Ci) emission in high-redshift (z>2) submillimeter galax- 
ies (SMGs) and quasar host galaxies (QSOs). Sensitive observations of the Ci(^Pi — >^ Po) and 
Ci(^P2 -^^ Pi) lines have been obtained at the IRAM Plateau de Bure interferometer and the IRAM 
30 m telescope. A total of 16 Ci lines have been targeted in 10 sources, leading to a total of 10 
detected lines — this doubles the number of Ci observations at high redshift to date. We include 
previously published C i observations (an additional 5 detected sources) in our analysis. Our main 
finding is that the C i properties of high-redshift systems do not differ significantly from what is 
found in low-redshift systems, including the Milky Way. The Ci(^P2 Pi)/Ci(^Pi ->-^ Po) and 
the Ci(^Pi -^^ Po)/"'^^CO(3-2) line luminosity (L') ratios change little in our sample, with respective 
ratios of 0.55±0.15 and 0.32dz0.13. The Ci lines are not an important contributor to cooling of the 
molecular gas (average -Lci/^fir ~ (7.7 ±4.6) x 10~^). We derive a mean carbon excitation temper- 
ature of 29.1±6.3K, broadly consistent with dust temperatures derived for high-redshift starforming 
systems, but lower than gas temperatures typically derived for starbursts in the local universe. The 
carbon abundance of X[C i]/X[H2]~(8.4±3.5) x 10~^ is of the same order as found in the Milky Way 
and nearby galaxies. This implies that the high-z galaxies studied here are significantly enriched in 
carbon on galactic scales, even though the look-back times are considerable (the average redshift of 
the sample sources corresponds to an age of the universe of ~2Gyr). 

Subject headings: galaxies: formation — cosmology: observations — infrared: galaxies — galaxies: 
starbursts — galaxies: evolution 



1. INTRODUCTION 

Detections of carbon monoxide (CO) and dust, in ob- 
jects at z >2 have enabled studies of the molecular gas 
properties in the early phases of galaxy formation (e.g., 
review by Solomon & Vanden Bout 2005). Molecular 
gas masses > 10"'^° Mq, enriched in C and O, provide 
strong evidence for powerful starbursts at early epochs 
in objects that may have later turned into massive spi- 
rals and/or elliptical galaxies. Until now, mid-J (J=3- 
6) CO lines and the continuum emission from dust have 
primarily been used to characterize the physical proper- 
ties of these massive gas reservoirs. One key approach is 
to observe the different rotational (J) transitions of CO 
(so-called 'CO line spectral energy distribution', or 'CO 
line ladder') and use radiative transfer models to con- 
strain molecular gas masses, densities, and temperatures 
(see the compilation by Weifi et al. 2007, 2011, hereafter 
Wll). In such analyses, however, it is frequently found 
that the density and kinetic temperature of the molecular 
gas are degenerate, so that data from additional tracers 
of the gas phase would be highly desirable. 

Such an alternative tracer of the cold molecular gas 
phase is atomic carbon (Ci). Studies of atomic carbon 
in the local universe have been carried out in molecular 
clouds of the galactic disk, the galactic center, M82 and 
other nearby galaxies (e.g.. White et al. 1994; Stutzki 
et al. 1997; Gerin & PhiUips 1998; 2000; Ojha et al. 
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2001; Israel & Baas 2002; Schneider et al. 2003, Israel 
2005). These studies have shown that Ci is closely asso- 
ciated with the CO emission independent of environment. 
Since the critical density for the Ci(l-O) and ^^CO(J — 
1 — J> 0) lines are both ~ lO^cm"^, this finding sug- 
gests that the transitions arise from the same volume and 
share similar excitation temperatures (e.g. Ikeda et al. 
2002). For reference, the critical densities for the Ci(l- 
0) and Ci(2-1) lines are 1x10^ cm~^ and 3x10^ cm~^, 
respectively (Schoier et al. 2OO5)0- As pointed out by Pa- 
padopoulos et al. (2004) such a close connection between 
C I and CO was not expected from early theoretical work 
(Tielens & Hollenbach 1985a, 1985b), where the C i emis- 
sion was thought to emerge from a narrow Cii/Ci/CO 
transition zone in molecular clouds. 

Because the '^P fine-structure system of atomic carbon 
forms a simple three-level system, detection of both op- 
tically thin carbon lines, Ci (^Pi -^^ Po) (rest frequency: 
492.161 GHz, hereafter C i(l-O)) and C i (^Pa -^'^ Pi) (rest 
frequency: 809.344GHz, hereafter Ci(2-1)), enables one 
to derive the excitation temperature, neutral carbon col- 
umn density and mass independently of any other in- 
formation (e.g., Ojha et al. 2001, Weifi et al. 2003). A 
combination of this method (using C i) with the afore- 
mentioned CO SEDs is particularly powerful as it elim- 
inates some of the degeneracy frequently found in CO 
radiative transfer models. Therefore, assuming that C i 
and CO trace similar volumes on galactic scales, the car- 
bon excitation temperature allows us to pin down the 
CO excitation conditions, as the excitation temperature 
of Ci and CO(l-O) are expected to be similar. 

For details see; |http:/ /www.strw.leidenuniv.nl/^ moldata/datafiles/catom.d 
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Fig. 1.— C 1(1-0) channel maps of SMMJ02399-0136. Each 
channel is 300kms~^ wide and the central velocity is given 
in the top left corner. Contours are shown in steps of Icr 
(1(7=0.31 mjy beam~^) 



Previous observations of Ci at high redshift (see ref- 
erences given in Section 4 and the related Table) have 
shown that both the Ci(l 0) and Ci(2-1) lines arc typi- 
cally significantly fainter than the respective CO lines at 
similar frequencies (partly because they arc typically op- 
tically thin). In this paper we report on observations of 
atomic carbon towards 10 high-redshift sources (a total 
of 16 C I lines) that encompass both submillimctcr galax- 
ies (SMGs) as well as quasar host galaxies (QSOs). These 
observations roughly double the number of high z C I line 
detections that have been presented in the literature pre- 
viously (in a dozen papers). In Section 2 we present the 
observations that have been obtainted with the IRAM 
Plateau de Bure Interferometer (hereafter: PdBI) as well 
as the IRAM 30 m telescope. In Section 3 wc discuss the 
observations of the individual sources. We compare the 
C I properties (including sources from the literature) in 
Section 4. Section 5 summarizes our findings. Through- 
out this paper we use a A cosmology with Hq = 71 
kms-i Mpc-\ Oa = 0.73 and flrn = 0.27 (Spergel etal. 
2003, 2007). 

2. OBSERVATIONS 

2.1. Source Selection 

Given the expected faintncss of the C i emission line 
we have selected sources that are bright in CO emis- 
sion. We have mostly targeted quasars using the IRAM 
30 m telescope as their linewidths are typically small (~ 
400kms~^), to minimize susceptibility to low-frequency 
baseline instablities in the single dish observations. We 
have typically used the (intrinsically more stable) PdBI 
in the case of the SMGs that are faint and have signifi- 
cant linewidths 800 — 1000 km s"-'^, e.g. Grcvc ct al. 
2005). The SMGs have been selected to be in the red- 
shift range 2.03<z<2.81 so that both Ci transitions can 
be observed in the 2 and 1 mm bands of the PdBI and 
30 m telescope. As discussed in the individual sources 
below, some of the SMGs host an AGN, i.e. the two pop- 
ulations that have historically been referred to as SMGs 
and quasars do overlap. Most of the quasars in our sam- 
ple have higher redshifts. The source names, coordinates, 
redshifts, targeted C i transitions and telescopes used are 
given in Table 1 and results for individual sources are dis- 
cussed in Section 3. 



2.2. IRAM Plateau de Bure Observations 

The data were taken in the years 2008 and 2009 

(project IDs S046 & S0B3) using the upgraded PdBI with 
its dual polarization receivers and a bandwidth of 1 GHz 
(single polarization). At 2 mm, typically five 15m anten- 
nas were available in the compact D configuration, with 
spacings from 24 to 86m. These short spacings had av- 
erage rms phase errors of 15° at 2 mm and 30° at 1mm. 
System temperatures were 150K in the lower sideband at 
2 mm, and 250 to 400K in the lower sideband at 1.3mm. 
The spectral correlators covered ~1500kms~^ at 2mm 
and 1000 km s^^ at 1mm, with instrumental resolutions 
of 8 and 4kms~^, respectively. The primary amplitude 
calibrators were 3C454.3 (variable, but typically 25 Jy at 
2mm and 13 Jy at 1.3mm), and MWC349 (nonvarying, 
with 1.9 and 1.6 Jy at 2 and 1.3mm). The typical un- 
certainties in the flux scales and overall calibration are 
^10% at 2 mm and ^-^15-20% at 1.3mm and should be 
added to the formal fitting errors quoted in the remainder 
of this paper. The observing program monitored phases 
every 20 min. The data processing program used water- 
vapour monitoring receivers at 22 GHz on each antenna 
to correct amplitudes and phases at 2 mm and 1.3 mm 
for short-term changes in atmospheric water vapor. All 
visibilities are weighted by the integration time and the 
inverse square of the system temperature. Source maps 
were made with 'natural weighting' of the visibilities to 
maximise the sensitivity. The tuned frequencies, the re- 
sulting beam sizes (and position angles), and the la rms 
over the full bandwidth (1 GHz) for each source are given 
in Table 2. The primary beam size of PdBI at 2 mm 
(1mm) is ~35" (~22"). Emission has been cleaned to 
the 1.50- level in a tight box around the source. In the 
case of non-detections no cleaning was performed. 

2.3. IRAM 30 m Observations 

Observations were carried out with the IRAM 30 m 
telescope in the years 2004-2007. We used the CD or 
AB receiver configurations with the C/D 150 receivers 
tuned to the C i(l-O) transition or the A/B 230 receivers 
tuned to C i(2-l). The beam size of the 30 m at 140 GHz 
(230 GHz) is f=s 17" (11"). Typical system temperatures 
were « 200 K and « 270 K during winter and summer re- 
spectively. The observations were carried out in wobbler 
switching mode, with a switching frequency of 0.5 Hz 
and a wobbler throw of 50" in azimuth. Pointing was 
checked frequently and was found to be stable within 
3". Calibration was obtained every 12 min using stan- 
dard hot/cold-load absorber measurements. Both Mars 
and Uranus were used for calibration. The antenna gain 
was found to be consistent with the standard value of 6.6 
(10) JyK-i at 140 (230) GHz (temperatures in units of 
T^). We estimate the flux density scale to be accurate 
to about 15%. 

Data were recorded using the 4 MHz filter banks on 
each receiver (512 channels, IGHz bandwidth, 4 MHz 
channel spacing) and were processed using the CLASS 
software. After dropping bad data only linear baselines 
were subtracted from individual spectra. The tuned fre- 
quencies and the noise in each 150kms~^ channel are 
given in Table 3. 
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Fig. 2.— Results for SMM J02399-0136. Left: Integrated Cl(l-O) einissiou from -500kins"i to +500kms"\ with contours starting at 
±2(T in steps of 2(t (1 (t=0.17 mjy). Middle: Continuum measurement at 129.2 GHz (line— free ehanncls of the C l(l-O) observations, using 
the same total bandwidth as for the line), same contouring and la sensitivity as in the left panel. No continuum is detected at 129.2 GHz. 
Right: Continuum emission at 212.5 GHz, with contours starting at ±2(j in steps of 2(7 (l<T=0.18mJy). The beam sizes are shown in the 
bottom left corner, respectively (see Table 2). 



TABLE 1 

Properties of Observed Targets 



source 


RA 


DEC 


z 


FWHM 


transition 


telescope 


references'' 




J2000.0 


J2000.0 




kms~^ 








SMM J02399-0136 


02:38:51.87 


-01:35:58.8 


2.808 


~1100 


Cl(l-O), CI(2-1) 


PdBI 


1, 2 


APM 08279+5255 


08:31:41.64 


+52:45:17.5 


3.911 


480+35 


CI(2-1) 


30 m 


3 


RXJ091 1+0551 


09:11:27.50 


+05:50:52.0 


2.796 


150+20 


Cl(l-O), CI(2-1) 


30 m 


4, 5 


SMM J123549+6215'" 


12:35:49.42 


+62:15:36.9 


2.202 


600+50 


Cl(l-O), CI(2-1) 


PdBI 


6 


BRI 1335-0417 


13:38:03.42 


-04:32:34.1 


4.407 


420+60 


CI(2-1) 


30 m 


7 


SMM J14011+0252 


14:01:04.93 


+02:52:24.1 


2.565 


190+11 


CI(2-1) 


30 m 


8, 9 


SMM J16359+6612 


16:35:54.10 


+66:12:23.8 


2.517 


500+100 


Cl(l-O), CI(2-1) 


30 m 


10, 11 


SMM J163650+4057'^ 


16:36:50.41 


+40:57:34.3 


2.385 


710+50 


Cl(l-O), C 1(2-1) 


PdBI 


12, 6 


SMM J163658+4105'^ 


16:36:58.17 


+41:05:23.3 


2.452 


800+50 


Cl(l-O), C 1(2-1) 


PdBI 


13, 6 


PSS J2322+1944 


23:22:07.25 


+19:44:22.1 


4.120 


200+70 


C 1(2-1) 


30 m 


14 


Note. — The coordinates and redshifts refer to the pointing center and redshift used for tuning the receivers. 



These values are taken from the references given in the last column. 
^ references (coordinates, redshifts, FWHM's) are; [1]: Fraycr et al. 1998; [2]: Genzel ctal. 2003; [3]: Downes ctal. 
1999; [4]: Hainline etal. 2004; [5]: Wll; [6]: Tacconi etal. 2006; [7]: Guilloteau et al. 1999; [8]: Frayer etal. 1999; 
[9]: Downes & Solomon 2003; [10]: Sheth etal. 2004; [11]: Kneib etal. 2005, [12]: Neri etal. 2003; [13]: Greve etal. 
2005; [14]: Cox etal. 2002. 
^ alternative name: HDF 76. 
alternative name: Elias N2 850.4. 
alternative name; Elias N2 850.2. 

(Figs. 7-12). 

3.1. PdBI Observations 
3.1.1. Cl Line Emission 

The main results and data analysis are described in 
the following: C i(l-O) emission (observed at 2 mm wave- 
length) has been detected at a significant level in three of 
the four sources that have been targeted. However, C i(2- 
1) emission (observed at 1mm wavelengths) has been 
tentatively detected only in one source (SMMJ02399- 
0136). The latter measurement is complicated by the 
fact that the faint Ci(2-1) emission is situated on top of 
bright continuum emission at 1 mm wavelengths. An ad- 
ditional complication is that the rest frequency of the 
CO(7-6) line (806.651GHz) is close to that of Ci(2- 
1) (809.344GHz), i.e. the velocity difference at Z--2.5 is 
(Al/)/!/ X c~1000kms~^ (where c is the speed of light). 
This means that the lines begin to overlap in frequency 
space at this redshift if the intrinsic full width zero in- 
tensity (FWZl) is larger than 1000 km s~^. This is not 
the case in most of our sources, but potentially compli- 
cates the determination of the continuum redward of the 
Ci(2-1) line. In the light of this we present conservative 



Fig. 3.— Cl(l-O) channel maps of SMM J123549+6215. Each 
channel is 300 km s^^ wide and the central velocity is given 
in the top left corner. Contours are shown in steps of Icr 
(1(7=0.43 mjybeam-l) 



3. RESULTS 

In the following we present the results for each source. 
The PdBI observations are summarized in Section 3.1 
(Figs. 1-6) and the 30 m observations in Section 3.2 
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Fig. 4.— Results for SMM J123549+6215. Le/t; Integrated Cl(l-O) emission from -400 km s'^ to +400 kms"!, with contours starting at 
±1(T in steps of Icr (1 cr=0.26 mjy). Middle: Continuum measurement at 153.3 GHz (line— free channels of the C 1(1-0) observations, using 
the same total bandwidth as for the line), same contouring and Itr sensitivity as in the left panel. No continuum is detected at 153.3 GHz. 
Right: Continuum emission at 252.8 GHz, with contours starting at ±2(7 in steps of 2cr (1 (t=0.20 mJy). The beam sizes are shown in the 
bottom left corner, respectively (see Table 2). 
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Fig. 5. — Results for SMM J163650+4057. Left: Integrated Cl(l-O) emission from -500kms~^ to +500 kms"'^, with contours starting 
at ±1(T in steps of Icr (1 (t=0.24 mJy). The Cl(l-O) line is not detected. Middle: Continuum measurement at 145.4GHz (line-free channels 
of the Cl(l— 0) observations, using the same total bandwidth as for the line), same contouring and Icr sensitivity as in the left panel. 
The continuum at 145.4 GHz is not detected. Right: Continuum emission at 239.1GHz, with contours starting at ±2(T in steps of 2tT 
(1 cr=0.22 mJy). The beam sizes are shown in the bottom left corner, respectively (see Table 2). 
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Fig. 6.— Results for SMM J163658+4105. Left: Integrated Cl(l-O) emission from -400 kms'^ to +200 kms'^, with contours starting 
at ±1.5cr in steps of 1.5cr (1 cr=0.29 mJy). Middle: Continuum measurement at 142.6GHz (line— free channels of the Cl(l-O) observations, 
using the same total bandwidth as for the line), same contouring and Icr sensitivity as in the left panel. There is a tentative detection of the 
142.6 GHz continuum in this source. Right: Continuum emission at 234.5 GHz, with contours starting at ±2cr in steps of 2cr (1 cr=0.22 mJy). 
The beam sizes are shown in the bottom left corner, respectively (see Table 2). 



5(7 upper limits for the C i(2-l) line fluxes. For the 2 mm 
continuum we give 3a upper limits. The a's used for de- 
riving the upper limits are based on the la noise in a 
channel that has the expected FWHM of the respective 
line. 

3.1.2. Continuum Emission 

Continuum emission at 1 mm has been detected in all 
cases (only one object is tentatively detected in the 2 mm 
band). We estimated the expected 1mm and 2mm con- 
tinuum flux densities (for the respective frequency) based 



on published SCUBA (850 /im) measurements of the in- 
dvidual sources (column 7 in Table 4), assuming a mod- 
ified black body (see more detailed discussion in Sec- 
tion 4.2.20- We use the best-fit values for SMGs by 
Kovacs et al. (2006, Tdust,SMG=34.6±3K, /3smg=1.5) 
and assign uncertainties of ±10 K to Tdust and ±0.2 to /?. 
In the case of the 1 mm observations, the measured flux 
densities of all our targets are consistent with these ex- 

^ We derive almost identical estimates using an Arp 220 spectral 
energy distribution 
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TABLE 2 

Plateau de Bure Observations 



source FWZP C I trans. i^ohs beamsize, P.A. rms'' Ici"^ Scont'^ Scxp"^ 

[kms-l] [GHz] "x",° [mjyb-i] [Jykms-l] [mjy] [mjy] 



SMM 02399-0136 


~1000 


Cl(l-O) 


129.2439 


4.9x3.5, 161 


0.11 


1.9+0.2 


<0.5 


1.2+0.5 






CI(2-1) 


212.5373 


3.0x2.0, 165 


0.18 


1.8+0.9 


6.2+0.2' 


7+3 


SMM 123549+6215 


~800 


C 1(1-0) 


153.7041 


4.6x3.3, 139 


0.17 


1.1+0.2 


<0.8 


0.8+0.3 






CI(2-1) 


252.762 


2.2x1.5, 145 


0.20 


<1.0 


1.7+0.2 


3.4+1.3 


SMM 163650+4057 


~1400 


Cl(l-O) 


145.3946 


4.1x2.8, 69 


0.17 


< 0.7* 


<0.8 


0.7+0.3 






CI(2-1) 


239.0966 


1.9x1.8, 99 


0.18 


<2.0 


3.7+0. 2S 


2.9+1.2 


SMM 163658+4105 


~1000 


Cl(l-O) 


142.5726 


3.7x2.8, 69 


0.16 


0.9+0.2 


~0.6'> 


0.8+0.3 






CI(2-1) 


234.456 


3.6x1.4, 174 


0.21 


<2.0 


2.0+0.2 


3.3+0.6 



Note. — The uncertainties given here are only based on the formal fitting error and do not include the uncertainties 
in flux calibration (as discussed in Sec. 2.2). 

The full width zero intensity (FWZI) has been estimated from published CO spectra (see references in Table 1). 
^ Rms in mjy beam" over the full 1 GHz bandwidth. 
Upper limits to the Cl(2— 1) line fluxes are 5cr values (see Section 3). 

Measured continuum flux density at frequency given in column 4. In the case of non-detections we quote 3<j upper 
limits. 

" Expected continuum at observed frequency based on SCUBA 850 fim measurements (Table 4, Column 7) and assuming 
a modifled blackbody with /3=1.5 and a dust temperature of T=34.6 K (Kovacs et al. 2006). The uncertainties reflect the 
model uncertainties of 40%, as discussed in Section 4.2.2. 

* Emission is marginally detected in the blue wing of the CO (3-2) line (Neri et al. 2003). 

s This source is potentially resolved in the continuum. The peak flux density is 2.8+0.2 mjy. 

^ This source is is only marginally detected in the 2 mm continuum (<3(t). 



TABLE 3 
IRAM 30 m Observations 



source 


C I trans. 


I'obs 


rms'' 


Fpcak 


FWHM 


Ici 






[GHz] 


[mjy beam"*] 


[mJy] 


[km s~*] 


[Jykms *] 


APM 08279+5255 


CI(2-1) 


164.7988 


1.6 (85kms) 


<4.8 




<1.1^ 


RXJ0911+0551 


C 1(1-0) 


129.6524 


2.9 (40kms) 


13.3 


140 + 25 


2.1+0.3 




CI(2-1) 


213.2096 


3.8 (45 kms) 


14.8 


145 + 30 


2.3+0.4 


BRI 1335-0417 


Cl(l-O) 


91.0229 


4.8 (60 kms) 


<14.6 




<2.2*' 




CI(2-1) 


149.6841 


1.7 (60 kms) 


<5.1 




<0.8'' 


SMM J14011+0252 


CI(2-1) 


227.0121 


2.3 (50 kms) 


14.4 


205 + 25 


3.1+0.3 


SMM J16359+6612 


CI (1-0) 


139.9229 


0.8 (75 kms) 


3.9,3.8 


250, leo'^ 


1.7+0.3 




C 1(2-1) 


230.0991 


0.9 (80 kms) 


5.4,4.4 


250, leo'^ 


2.2+0.3 


PSS J2322+1944 


C 1(2-1) 


158.0780 


3.2 (60 kms) 


8.0 


160 + 55 


1.4+0.3 


Velocity resolution 


for given rms 


shown in 


brackets. 









^ 3(7 upper limit over the full width of the expected line width. 
Profile fit with two Gaussian fits (e.g. Weifi et al. 2005a). 



trapolations and inspection of individual datacubes did 
not reveal the presence of an additional C i(2-l) line. No 
significant continuum emission was detected in the 2 mm 
data, consistent with expectations (Table 2). 

All derived line fluxes and continuum flux densities are 
summarized in Table 2. In the following we briefly dis- 
cuss the individual sources. 

3.1.3. SMM J02399-0136 

This source is one of the first SlVIGs 
(S850/im=26±3mJy, Ivison et al. 1998) ever detected in 
CO emission (J=3, Frayer et al. 19980. It has subse- 
quently been studied in detail in CO emission (J=3: 
Genzel et al. 2003, J=l: Ivison et al. 2010a). Ivison et 
al. (1998) noted that this source hosts a dusty AGN; its 
QSO lines are discussed in Villar-lVIartin et al. (1999) 
and Vernet & Cimatti (2001). In the LpiR-ix-ray plane 
S1VI1VIJ02399-0136 has the same luminosities as the 
Cloverleaf (Alexander et al. 2005a). Ivison et al. (2010b) 

^ J=3: CO(3-2) 



conclude that this source comprises a merger between a 
FIR-luminous starburst, a QSO host and a faint third 
component. 

For our analysis we adopt a FWZI of ^ 1000 km s^^ for 
SMM J02399-0136. The Ci(l-O) line is clearly detected 
in both individual channel maps (Figure 1) and the in- 
tegrated line emission (Figure 2, left). The total flux of 
the Ci(l-O) hue is 1.9±0.2 Jy kms^^. The continuum at 
2 mm is not detected at a 3(7 limit of Si29ghz <0.51mJy 
(Figure 2, middle). The source is clearly resolved in the 
1 mm continuum at 212.5 GHz with a total flux density of 
S2i2GHz=6.2±0.2mJy (Figure 2, right). Given the large 
hnewidths in SMM j62399-0136 it is difficuh to separate 
the continuum from possible Ci(2-1) emission in this 
source at 1 mm wavelengths. Genzel et al. (2003) have 
derived a total continuum flux of 7.0±1.2mJy measured 
at 235 GHz. If we assume that this is the correct contin- 
uum flux, we would expect a 212 GHz flux of ~4.9mJy, 
or about 1.3 mJy less than what is observed. We con- 
clude that some of the measured flux may be attributed 
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to C i(2-l) emission in this source and derive a tentative 
Ci(2-1) flux of '-^l.S Jykms"^ with a considerable error 
bar (50%). We note that SMM J02399-0136 is thought 
to be lensed (magnification factor fj,—2.5, Frayer et al. 
1998, Genzel et al. 2003). 

3.1.4. SMMJ123549+6215 

This source is sometimes also referred to as HDF 76 
(S850/xm=8.3±2.5mJy, Chapman et al. 2003) and has 
been detected in CO emission (J=3 and J— 6: Tac- 
coni et al. 2006, 2008, J=l: Ivison et al. 2010a). 
SMM J123549+6215 is X-ray bright (Alexander et al. 
2005a) and its luminosity is classified as AGN-dominated 
(see also Alexander et al. 2005b, Takata et al. 2006). 

We adopt a FWZI of -800 km s^^ for 
SMMJ123549+6215. The Ci(l"0) line is detected 
both in individual channel maps (Figure 3) and in the 
integrated line emission (Figure 4, left) with a total flux 
of l.l±0.21Jy kms~^. The continuum at 2mm is not 
detected at a 3a flux density limit of Sissghz <0.78mJy 
(Figure 4, middle). The 1 mm continuum emission is de- 
tected at S252.8GHz=l-7±0.2 mJy (Figure 4, right). This 
is only marginally consistent with the tentative detection 
quoted by Tacconi et al. 2006 (S213 GHz=2.0±0.6mJy, 
the extrapolated flux at 213 GHz would be —1 mJy), and 
within the expectations based on the extrapolation from 
the SCUBA measurement (Table 2). We derive a 5a 
upper limit of the Ci(2-1) line (assuming a bandwidth 
of 800kms"i) of <1.0Jykms-^ 

3.1.5. SMMJ163650+4057 

This source is also known as Elias N2 850.4 
(S850pm=8.2±1.7mJy, Ivison et al. 2002) and first CO 
data {J— 3 and J=7) have been presented by Neri et al. 
(20030 ( J=l data are shown in Ivison et al. 2010a). Im- 
proved CO data are shown in Tacconi et al. (2006, 2008, 
J— 3 and J— 6). The AGN broad lines of this source 
(which in fact is a pair of blue and red galaxies, Ivison et 
al. 2002) are discussed in Swinbank et al. (2005, see also 
Takata et al. 2006, Menendez-Delmestre 2007, 2009). 

The Gi(l-O) line is not detected in our data 
with a 3a upper limit across the expected linewidth 
(~1000kms-i) of 0.7Jykms-i (Figure 5, left, but 
see caption of Table 2). Likewise, the continuum at 
2 mm is not detected at a 3a flux density limit of 
Si45GHz <0.74mJy (Figure 5, middle). The 1mm con- 
tinuum emission appears to be slightly extended with 
S239.iGHz=3.7±0.4 mJy (peak flux: 2.8±0.18mJy, Fig- 
ure 5, right), in rough agreement with our SCUBA ex- 
trapolation and with Neri et al. (2003) and Tacconi 
et al. (2006) who quote consistent 1mm flux densities 
for this source (S230 GHz=2.5±0.4mJy and 2.6±0.5mJy, 
respectively). We derive a 5a upper limit of the 
Ci(2-1) line (assuming a bandwidth of 1000 kms~^) of 
<2.0 Jykms^i. 

3.1.6. SMMJ163658+4105 

This source is also known as Elias N2 850.2 
(S850Mm=10-7±2.0mJy, Ivison et al. 2002) and CO data 
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Fig. 7.— Left: C 1(2-1) spectrum of APM 08279+5255; the line 
is not detected. 

are published in Greve et al. (2005, J=3)0, Tacconi et al. 
(2006, J=3 and J=7) and Ivison et al. 2010a (J=l). 

The Ci(l-O) line is detected at 0.90±0.18 Jy kms'^ 
over a velocity width of 600 kms~^ (Figure 6, left). This 
is the only source in which the 2 mm continuum emis- 
sion may be marginally detected at Si42.6GHz ~0.6mJy. 
If real (as suggested by the good match to the extrap- 
olated SCUBA flux density. Table 2), the continuum 
emission would contribute —40% to the Ci(l-O) flux 
density. Given the low S/N of the potential contin- 
uum flux, we do not apply this correction to the Ci(l- 
0) flux in the analysis that follows but note that such 
a correction would not change the statistical analysis 
below. The 1 mm continuum emission is detected at 
S234GHz=2.0±0.22 mJy (Figure 6, right), in good agree- 
ment with the extrapolation from the SCUBA measure- 
ment (and consistent with the estimate of Tacconi et al. 
2006, of S233GHz=l-5±0.5mJy, their Table 1). 

3.2. 30 m observations 

We have observed both the C i(l-O) and C i(2-l) emis- 
sion in 6 additional sources, four quasars and two sub- 
millimeter galaxies. All fitted parameters and line fluxes 
are summarized in Table 3. In the following we briefly 
discuss the individual sources. 

3.2.1. APM08279+5255 

This source is one of the (apparently) brightest 
sources known to date, with a submillimeter flux of 
S850Mm=75±4mJy (Lewis et al. 1998). The first CO 
detections by Downes et al. (1999, J=4, J=9) already 
indicated the extreme excitation conditions of the molec- 
ular gas (studied in detail in Weifi et al. 2007). One of 
the explanations for its apparent brightness and extreme 
excitation is that the lensing geometry is such that the 
central 100 pc of this source is magnified by a large fac- 
tor (/i=60-100 Weifi et al. 2007, see also Egami et al. 
2OOOB The Ci(l-O) line has been detected by Wagg et 

* In Greve et al. (2005) the source is referred to as 
SMM J16366-I-4105 

^ We note that Riechers et al. 2009a derived a magnification 
factor of only ^^4 for the molecular gas phase. 
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Left: C 1(1-0) spectrum of RX J0911+0551. Right: C 1(2-1) spectrum of RX J0911+0551. See Table 3 for respective fluxes and 
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Fig. 9. — Cl upper limits for BRI 1335-0417. Left: C 1(1-0) observations. Right: C 1(2-1) observations including the adjacent CO(7-6) 
line. The arrows in both panels indicate the expected central position of the C I lines. 



al. (2006). The source is not detected in Ci(2 1) emis- 
sion with an upper Hmit of IJykms"^ (Figure 7, left; 
upper limits are given in Table 3). 

3.2.2. RXJ0911+0551 

This source is bright in the submillimeter 
(S850/im=26.7±1.4mJy, Barvainis & Ivison 2002) 
and CO has been tentatively detected in this source by 
Hainline ct al. (2004, J=3). Follow-up observations by 
Wll (J=3,5,7,8,9) show that the CO emission line is 
very narrow (FWHM~110kms~^). We detect both the 
Ci(l-O) and Ci(2-1) emission lines in this source (see 
Figure 8 and Table 3). 

3.2.3. BRI1335-0417 

BRI1335-0417 is bright in the submillimeter (e.g., 
Benford et al. 1999, we adopt S850/^m=22±4mJy from 
their Figure 2) and has first been detected in CO ( J=5) 
by Guilloteau et al. (1997). High angular resolution 



CO observations revealed a complex, possibly merg- 
ing system (Carilli et al. 2002a, Riechers et al. 2008a, 
J=2). The source is not detected in both the Ci(l-O) 
and Ci(2 1) line with upper limits of 2.2 Jykms"^ and 
0.8Jykms"i (Figure 9, right; see also Table 3). The 
spectrum of the Ci(2-1) Unc in Figure 9 also includes 
the CO (7-6) line that is close in frequency space (the 
CO emission is discussed in Wll). 

3.2.4. SMMJ14011+0252 

SMMJ14011-I-0252 is one of the first submillimeter 
galaxies (S850/im=15±2mJy, Ivison et al. 2000) that 
has been detected in CO emission (Frayer et al. 1999, 
J=3; Downes & Solomon 2003, J=3, J=7). Downes & 
Solomon argue that the source is strongly lensed, with a 
magnification factor /i=25±5. The C i(2-l) line is shown 
in Figure 10 (see also Table 3). 

3.2.5. SMMJ16359+6612 
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Fig. 10.— C 1(2-1) spectrum of SMM J14011+0252. See Ta- 
ble 3 for fluxes and noise (see WeiB et al., 2005b, for tfie C 1(1-0) 
spectrum of the same source) . 

This source is a triple-lensed system in which ah three 
images are bright in the submilhmeter (brightest image 
'B': S850/.m = 17±2mJy, Kneib et aL 2004). All images 
have subsequently been detected in CO emission (Sheth 
et al. 2004, J=3, Kneib et al. 2005, J=3). Kncib et al. 
(2004) have argued that the magnifications for the three 
images are 14±2, 22±2 and 9±2, respectively. Higher- 
J CO transitions have been presented in Weil3 et al. 
(2005a). The beams are 18" and 11" at our tuned fre- 
quencies at 2 mm and 1 mm respectively, i.e. only the 
brighest image ('B') is covered in our observations. In 
Figure 11 we show the Ci(l-O) (left) and the Ci(2~l) 
(right) spectra (see also Table 3). 

3.2.6. PSSJ2322+19U 

PSSJ2322+1944 is a lensed source (//=5.4±0.3, 
Riechers et al. 2008b) which is bright in the sub- 
millimeter (S850/im=22.5±2.5mJy, Isaak et al. 2002; 
Si3ooMm=9.6±0.5mJy, Omont et al. 2001). Strong CO 
emission has been reported by Cox et al. (2002, J=4, 
J—b) and resolved imaging has revealed a molecular Ein- 
stein ring in this source by Carilli et al. (2002b, 2003, 
J=2) and Riechers et al. (2008b, J=2). The Ci(l-O) 
line has been detected by Pety et al. (2004). We detect 
the Ci(2-1) line at high significance (Table 3, Figure 12 
shows a wider bandwidth spectrum after combining the 
C i(2-l) data with the CO(7-6) data presented in WU). 

4. DISCUSSION 

4.1. Previous high-z Ci observations 

We have observed C i emission in 10 high-redshift sys- 
tems, targeting a total of 16 Ci lines, and have detected a 
total of 10 lines (one more tentatively). In the discussion 
that follows we include C i measurements of high-redshift 
sources from the literature (an additional 7 sources) . All 
C I measurements (including upper limits) to date are 
summarized in Table 4. We have re-calculated all line 
and far-infrared luminosities for the published sources 
to ensure that consistent assumptions and cosmological 
parameters are used in the analysis that follows. 



4.2. C I, CO and FIR luminosities 
4.2.1. Line luminosities 

In the following we derive line luminosities both in 
units of Lq (Liinc) and Kkms~^ pc^ (^line)' where 'line' 
refers to either CO(3-2), Ci(l-O) or Ci(2-1). The for- 
mer luminosity is given by 

Lnnc^lMxlO-^SiinoAvi^,,,t{l + z)-^Dl, (1) 

(Solomon et al. 1992), where the luminosity of the 
line, Liino, is measured in Lq; the velocity integrated 
flux, S'linoAu, in Jykms"^; the rest frequency, i^rcst = 
t'obs(l + z), in GHz; and the luminosity distance, Z^l, in 
Mpc (we have tabulated the luminosity distances in Ta- 
ble 5 assuming the cosmological parameters given at the 
end of Section 1). Luno (Eq- 1) is a true luminosity, i.e. 
it presents the amount of cooling that is radiated away 
by the respective line. 

The line luminosity can also be expressed in units of 
Kkms^^ pc^: 

LL, = 3.25 X 10^ Sunc >^-l Dl (1 + z)-\ (2) 

(Solomon et al. 1992). As noted by Solomon & Vanden 
Bout (2005), as Lfj^^j, is proportional to the brightness 
temperature, the L' ratio of two lines is a measure of the 
ratio of their intrinsic brightness temperatures (assuming 
that both lines arise from the same volume, i.e. same 
area and linewidth). 

We have calculated both luminosities for the Ci(l~0), 
Ci(2-1) and CO(3-2) lines and they are given in Table 5 
(uncorrected for possible magnification). 

4.2.2. Far-Infrared (FIR) luminosities 

To derive consistent results within our sample, we have 
recalculated FIR luminosites for all sources based on 
their measured 850/im flux (as tabulated in Table 4). 
The reason why we chose the 850/im band as our refer- 
ence wavelength is because continuum emission has been 
measured in this band for almost every source in our 
sample. 

Following Blain et al. (2003), a simple description of 
the far-infrared SED as a function of frequency v is based 
on a backbody spectrum Bi, ^ / [exp{hv / kT) — 1] at 
temperature T, modified by a frequency-dependent emis- 
sivity function ^ v^^ ('modified black body'). This 
yields a functional behavior at the source's rest frame of 
^ u^+P l[exp{hvlkT) ~ 1] (e.g. Blain et al. 2003). 

At wavelengths where the dust is optically thin, the 
observed flux is proportional to the dust mass and the 
the observed flux density can be written as: 

■^"o = ^''^^'2^'* A^dustK(l'r)B^,(Tdust) 

(Downes et al. 1992), where k{v) is the dust absorption 
coefficient and the subscripts o and r refer to the observed 
and restframe frequencies. 

We have derived FIR luminosities by integrating under 
a modified blackbody curve from 40 to 400 [im. (Sanders 
et al. 200313- The numbers are given in the last col- 
umn in Table 5. For the quasars in our sample we have 

The FIR luminosities decrease by 20-30% if the definition by 



TABLE 4 

Compilation of all high-z C i observations to date 
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z 


magn. (//) 
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Ici(l-O) 


lci{2-l) 




references for columns'' 














[Jy kms~^] 


[Jy kms~'] 


[Jy kms~'] 


[mJy] 






U] 


[51 


[6] 


[7] 


[1] 


[2] 


[3] 


[4] 


[5] 


[6] 


p] 














SMM J02399-0136 


2.808 


2.5 


3.1+0.4 


1.9+0.2 


1.8+0.9 


26+3 


1,2 


1,2 


2 
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APM 08279+5255 


3.911 


60-100'' 


2.2+0.2<: 


0.93+0.13 


<1.1 


75+4=1 
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5 
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7,8 


RXJ0911+0551 


2.796 


~20 


2.0+0.3= 


2.1+0.3 


2.3+0.4 


26.7+1.4 


9,10 


8 


9,10 
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F10214 


2.285 


IQf 


4.2+0.8 


2.0+0.4 


4.6+O.7S 


50+5 


11 


11 


11 


12,13 


13 
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SDSS J1148+5251 


6.419 


1 


0.18+0.04 




0.22+0.05 


7.8+0.7 


14,15 




15,16 




17 


18 


SMMJ123549+6215 


2.202 


1 


1.6+0.2 


1.1+0.2 


<1.0 


8.3+2.5 


19,20 




19,20 








21 


GN 20.2 


4.051 


1 


0.5+0.2== 




<1.9'^ 


~9.9+2.3 


22 


23 


22 




24 


25 


GN 20 


4.055 


1 


0.8+0.1== 




<1.2h 


20.3+2.1 


22 


23 


22 




24 


25 


BRI 1335-0417 


4.407 


1 


1.0+0.2' 


< 2.2 


<0.8 


23+3j 


26 
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27 


SMM J14011+0252 


2.565 


25+5!^ 


2.8+0.3 


1.8+0.3 


3.1+0.3 


15+2 


28,29 


29 


28,29 


12 





30 


Clovcrlcaf 


2.558 


11' 


13.2+0.2 


3.9+0.6 


5.2+0.3 


58.8+8.1 


12 


31 


32 


12.33 


32 


8 


SMM J16359+6612 


2.517™ 


22±2 


2.8+0.2 


1.7+0.3 


1.6+0.3 


17+2 


12,34,35 


35 


12,34.35 








36 


SMM J163650+4057 


2.385 


1 


1.6+0.2 


<0.72 


<2.0 


8.2+1.7 


37.19,20 




37.19,20 








38 


SMMJ163658+4105 


2.452 


1 


1.8+0.2 


0.90+0.2 


<2.0 


10.7+2.0 


39,19,20 




37,19,20 








38 


MM 18423+5938 


3.930 


20° 


2.8 + 0.3c 


2.3 + 0.5 


4.2 + 0.8 


70+5° 


40 




40 


40 


40 


40 


SMMJ213511-0102 


2.326 


32.5 + 4.5 


13.2 + 0.1 


16.0 + 0.5 


16.2 + 0.6 


113 + 13 


41 


42 


43 


43 


43 


42 


PSS J2322+1944 


4.120 


5.3+0.3 


2.1+0.3== 


0.80+0.12 


1.4+0.3 


22.5+2.5 


44 


45 


44 


46 





47 



References in last column are the following papers: [0] this work; [1] Prayer et al. 1998; [2] Gcnzel ct al. 2003: [3] Ivison ct al. 1998: [4] Downes ct al. 1999; 
[5] WeiB et al. 2007; [6] Wagg et al. 2006; [7] Lewis ct al. 1998; [8] Barvainis & Ivison 2002; [9] Wll; [10] Hainline et al. 2004; [11] Downcs et al. 1995; [12] 
WeiB et al. 2005b; [13] Ao et al. 2008; [14] Bertoldi et al. 2003; [15] Walter et al. 2003; [16] Walter et al. 2004; [17] Riechers et al. 2009b; [18] Robson et al. 
2004; [19] Tacconi et al. 2006; [20] Tacconi et al. 2008; [21] Chapman et al. 2003; [22] Daddi et al. 2009a; [23] Daddi et al. 2009b; [24] Casey et al. 2009; 
[25] Pope et al. 2006; [26] Guilloteau et al. 1997; [27] Benford et al. 1999; [28] Prayer et al. 1999; [29] Downes & Solomon 2003; [30] Ivison et al. 2000; [31] 
Venturini & Solomon 2003; [32] Weifi et al. 2003; [33] Barvainis et al. 1997; [34] Sheth et al. 2004; [35] Kneib et al. 2005; [36] Kneib et al. 2004; [37] Neri et 
al. 2003; [38] Ivison et al. 2002; [39] Greve et al. 2005; [40] Lestrade et al. 2010; [41] Swinbank et al. 2010; [42] Ivison et al. 2010c; [43] Danielson et al. 2010; 
[44] Cox et al. 2002; [45] Riechers et al. 2008b; [46] Pety et al. 2004; [47] Isaak et al. 2002. 

^ An average value of 80 is used for further analysis. See Riechers et al. (2009a) for a significantly lower fj, for the CO emission. 

= No CO(3-2) measurements available. CO(3-2) fluxes have been derived from measured CO(4-3) fluxes assuming that the emission is thermalized up to 

J=4 (i.e., Sco{3-2) = (|)^xSco(4-3)- 

=1 Barvainis & Ivison give S850Mm=84+3 mJy. 

= Data from Wll. Hainline et al. (2004) reported a significantly larger velocity width/total flux. 

' Downcs etal. (1995) derive slightly different magnifications for the CO (^=10) and FIR (;u=13) emission. 

s Consistent with the 2cr upper limit of 7 Jykms^' given by Papadopoulos (2005). 

2(7 upper limits by Casey et al. (2009). 
' Extrapolated from the J=5 measurement by Guilloteau et al. (1997). 
j Value taJcen from Figure 2 of Benford et al. (1999). 

^ Ivison etal. (2001) adopted a magnification value that is lower by an order of magnitude (/U=2.5). 
' Kneib etal. (1998) give a magnification range of /l(=18-30. 
™ Only the brightest component ('B') is covered here. 

° Magnification fa<;tor unknown (Lestrade et al. 2010). We here assign a fa<;tor of 20 to this source; consistent with the large fa<;tors found for other strongly 
lensed SMGs. 

° The 850 micron flux is derived from the 1.2 mm flux given in Lestrade et al. 2010, assuming a modifled black body with Tdust=40K and ,9=1.6. 
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Fig. 11.— Left: Cl(l-O) spectrum of SMM J16359+6612, Right: Cl(2-1) spectrum of SMM J16359+6612. Two Gaussian functions have 
been fitted to account for the two components (e.g., Sheth et al. 2004, Kneib et al. 2005, Weifi et al. 2005a). See Table 3 for fluxes and 
noise. 
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Fig. 12.— C 1(2-1) spectrum of PSS 2322+1944 with the adjacent 
CO(7-6) line (Wll). The velocity and redshift scales refer to the 
C 1(2-1) observations. See Table 3 for fluxes and noise. 

used the best-fit values (Tdust=47±3K, /3=1.6±0.1) de- 
rived by Beelen et al. (2006). Likewise, we have used 
the best-fit values by Kovacs et al. (2006) for the SMGs 
(Tdust=34.6±3K, P—l.b). In both cases, following Bee- 
len et al. (2006) and Kovacs et al. (2006) we assign 
uncertainties of ±3K to Tdust and ±0.1 to $. Based on 
this we assign uncertainties of A(log Lfir)=0.15, corre- 
sponding to ^ 40%, to our FIR luminosity estimates (the 
measurement uncertainty in the 850/xm fluxes is signifi- 
cantly lower than this number). 

All luminosities given in Table 5 have not been cor- 
rected for magnification. To derive the unlensed lumi- 
nosities, all numbers need to be divided by the respective 
magnification factors (/i) given in Table 4. 

4.3. Ci, CO and Ci/CO line ratios 

In Figure 13 we plot the C i(l-O) and C i(2-l) line lumi- 
nosities as a function of CO(3-2) luminosities. Most Ci 
luminosities are between 10% and 100% of the CO (3-2) 
luminosities (as indicated by the dashed diagonal lines). 
We calculate the -^'ci(i-o)/-^co(3-2) ratios of all targets 
(Table 6). We chose the CO(3-2) transition as a mea- 
sure for the CO luminosity as this transition has been 
observed for almost all targets. We find an average value 
of Lq^^_p-|/Lqq^3_2)=0.32±0.13 with no significant dif- 
ferences between the QSO and SMG sample. We plot 
the -^ci(i-o)/-^co(3-2) ratios in Figure 14 as a function 
of .^ci(i-o)/-^FiR (discussed in Section 4.7). 

We now compare the -^ci(i-.o)/-^co(3-2) ratios to sys- 
tems at low-redshift (including our Galaxy). This 
comparison is complicated by a two facts: (a) in 
most low-z cases, the CO luminosity has been derived 
from the J=l transition, not J=3 as in this study. 
We here assume -^^co(i-o) 
et al. 2005a, Ao et al. 2008) 



0-9 X L'co(3-2)' Weifi 
- this translates to a 



-^Ci(i-o)/-^CO(i-o)~^-^9±0.12 for the high-redshift sam- 
ple; (b) The beam sizes of the Ci(l-O) and CO obser- 
vations are typically different for the low-z sources (this 
fact does not play a role at high redshift, where our beam 
sizes are in most cases larger than the source size). 



Gerin & Phillips (2000) found an average value of 
-^ci(i-o)/-^co(i-o)~'-'-2^'-'-^ iri nearby galaxies and sim- 
ilar values have been derived in other studies: e.g., M83: 
0.18±0.04, NGC6946: 0.20±0.04 (Israel & Baas 2001, 
their table 4), IC342: 0.16±0.03, Maffci 2: 0.10±0.05 
(Israel & Baas 2003, their table 4), M51: 0.22±0.06 (Is- 
rael et al. 2006, their table 4), Hennize2-10: 0.09±0.03, 
NGC 253: 0.22±0.05 (Bayet et al. 2004, their tables 2 
and 3), M82: 0.^3- The average of the same value in 
the Milky Way is around 0.15±0.1; Fixsen et al. (1999) 
(0.09, Ojhaet al. (2001): 0.08, Okaet al. (2005): ranging 
from 0.05-0.3). 



We conclude that, to first order, the L[ 



CI(l-0)/-^CO 

tio does not appear to vary strongly with either envi- 
ronment or redshift, from low-excitation environments 
like the Milky Way all the way to extreme starbursts in 
the early universe. This finding is in agreement with our 
earlier conclusion based on a significantly smaller sample 
(Weii3 et al. 2005b). 

4.4. Atomic Carbon Excitation Temperature 

Some sources have now been detected in both C i lines 
and we can use their flux ratio to derive the carbon ex- 
citation temperature. As discussed in Schneider et al. 
(2003) the carbon excitation temperature in the local 
thermodynamic equilibrium (LTE) can be derived from 
the C I line ratio via the formula 



Te, = 38.8 X ln{—)-\ 
K 



(3) 



where ^=ici(2-i)/^ci(i-o)- This assumes that both 
carbon lines share the same excitation temperature and 
are optically thin. The latter assumption can be tested 
using the equations given in Schneider et al. (2003, 
their equations A. 6 and A. 7). Since these equations 
require knowledge of Tox we assume here, as a first 
guess, Tcx=Tdust (Sec. 4.2.2) and intrinsic source sizes 
of r=lkpc. This yields optical depths for the Ci(l-O) 
hue between 0.05-0.3, justifying the use of equation 3. 
Larger source sizes, as discussed by Tacconi et al. (2008), 
will lead to even lower optical depths. 

We summarize the derived excitation temperatures in 
Table 6 and find values between ^25 and '^45K with 
an average value of 29.1±6.3K. Note that if we used 
these derived Tox, the C i lines are still optically thin. In 
the local universe, the carbon excitation tempertature 
could only be derived for a few objects to date (using 
the same method as employed here): M82: ^55 K (us- 
ing the values given in Stutzki et al. 1997), NGC 253: 
-20 K (Bayet et al. 2004), Milky Way: 17.5 K (Fixsen et 
al. 1999). Herschel SPIRE FTS observations wih soon 
provide measurements of both C i lines in a sample of 
nearby ULIRGS (Van der Werf et al. 2010). 

4.5. Atomic Carbon Masses 

We now derive the atomic carbon masses of our targets. 
As derived in Weif5 et al. (2003, 2005b), the carbon mass 
can be calculated from the Ci(l-O) line luminosity via: 

Using the Cl(l-O) measurement from Stutzki et al. 1997 and 
the CO(l-O) data presented in Walter et al. 2002. 
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TABLE 5 

Derived Luminosities based on the values given in Table 4. 



source 




T b 

iJCI(l-O) 


Lci(2-l)'' 


Lco(3-2)'' 


'^CI(l-O)'' 


-'^01(2-1) 


'^00(3-2)'' 


Lfir'' 




[Gpc] 


[10^ Lq] 


[lO^L©] 


[IO^Lq] 


[10i"Kkms-ipc2] 


[lO^OKkms- 


lpc2] [lO^^Kkms-ipc^] 


10^3 Lo 


SMMJ02399-0136 


23.84 


14.5±1.5 


22±11 


17±2 


3.8±0.4 


1.3±0.7 


12.6±1.6 


1.41 


APM08279+5255 


35.55 


12.3±1.7 


<24 


20.4±1.9 


3.2±0.5 


<1.4 


15.4±1.4 


9.45 


RXJ0911+0551 


23.72 


16±2 


29±4 


10.7±1.4 


4.2±0.6 


1.7±0.3 


8.1±1.1 


3.81 


F10214 


18.53 


10.9±1.9 


40±6 


16±3 


2.9±0.5 


2.4±0.4 


12±2 


7.58 


SDSSJ1148+5251 


63.79 




10±2 


3.6±0.8 




0.60±0.14 


2.7±0.6 


0.87 




17.71 


5.5ibl.l 


<8.2 


5.6i0.7 


1.5i0.3 


<0.5 


4 3±0.5 


0.45 


GN20.2 


37.08 




<44 


5±2 




<2.6 


3.7±1.5 


0.55 


GN20 


37.12 




<28 


7.8±1.0 




<1.6 


5.9±0.7 


1.12 


BRI1335-0417 


40.99 


<35 


<21 


11±2 


<9.2 


<1.2 


8.5±1.7 


2.78 


SMMJ14011+0252 


21.35 


12±2 


33±3 


12.9±1.4 


3.1±0.5 


1.97±0.19 


9.7±1.0 


0.82 


Cloverleaf 


21.28 


25±4 


56±3 


60.4±0.9 


6.7±1.0 


3.28±0.19 


45.7±0.7 


8.64 


SMMJ16359+6612 


20.86 


10.8±1.9 


23±3 


12.5±0.9 


2.8±0.5 


1.35±0.18 


9.4±0.7 


0.93 


SMMJ163650+4057 


19.53 


<4.2 


<19 


6.5±0.8 


<1.1 


<1.1 


4.9±0.6 


0.45 


SMMJ163658+4105 


20.21 


5.5±1.1 


<20 


7.7±0.9 


1.4±0.3 


<1.2 


5.8±0.6 


0.58 


MM18423+5938 


35.76 


30 ±7 


92± 17 


26 ±3 


8.0 ± 1.7 


5.4 ± 1.0 


20 ±2 


3.84 


SMMJ213511-0102 


18.94 


88 ± 3 


147 ±5 


51.2 ±0.4 


23.2 ±0.7 


8.7 ±0.3 


38.7 ±0.3 


6.15 


PSSJ2322+1944 


37.83 


11.5±1.7 


33±7 


21±3 


3.0±0.5 


1.9±0.4 


16±2 


2.78 


Note. — All luminosities 


given in this table are uncorrected for lensing. To correct for gravitational 


magnification, all numbers have to be 



divided by fj, given in Table 4, Column 3. 

For Ho = 71 kms-i Mpc-\nA = 0.73 and = 0.27 (Spergel etal. 2003, 2007). 
^ From Equation 1, all units in lO'^ Lq. 

From Equation 2, all units in 10^" Kkms^^ pc^. 

See discussion in Section 4.2.2. 



TABLE 6 
Derived Physical Properties. 



source 


a 

-1- cx 


Mci 


Mh2 


X[C I]/X[H2] 


Lci(i-o)/Lfir 


'^CI(l-0)/'^'cO(3-2) 


T ' /T ' 

CI(2-l)/^CI(l-0) 




[K] 


[10^ Mq] 


[10"' Mq] 


10-5 


10-6 


SMMJ02399-0136 


21.6±6.3 


2.0±0.2 


4.0±0.5 


8.1±1.4 


10±4 


0.30±0.05 


0.35±0.18 


APM08279+5255 


30.00 


0.050±0.007 


0.154±0.014 


5.4±0.9 


1.3±0.6 


0.21±0.03 


<0.44 


RXJ0911+0551 


23.5±3.0 


0.26±0.04 


0.32±0.04 


14±3 


4.2±1.8 


0.52±0.10 


0.41±0.08 


F10214 


42.0±10.9 


0.36±0.06 


0.95±0.18 


6.4±1.7 


1.4±0.6 


0.24±0.06 


0.8±0.2 


SDSSJ1148+5251 


30.00 




2.2±0.5 










SMMJ123549+6215 


30.00 


1.8±0.3 


3.4±0.4 


9±2 


12±5 


0.34±0.08 


<0.34 


GN20.2 


30.00 




3.0±1.2 










GN20 


30.00 




4.7±0.6 










BRI1335-0417 


30.00 


<11 


6.8±1.4 


<28 


<13 


<1.1 


<0.13 


SMMJ14011+0252 


32.4±5.2 


0.15±0.03 


0.31±0.03 


8.2±1.6 


14±6 


0.32±0.06 


0.64±0.12 


Cloverleaf 


26.7±3.0 


0.76±0.12 


3.32±0.05 


3.8±0.6 


2.9±1.3 


0.15±0.02 


0.49±0.08 


SMMJ16359+6612 


26.2±4.0 


0.16±0.03 


0.34±0.02 


7.8±1.5 


12±5 


0.30±0.06 


0.48±0.11 


SMMJ163650+4057 


30.00 


<1.4 


3.9±0.5 


<5.8 


<9 


<0.22 


<1.0 


SMMJ163658+4105 


30.00 


1.8±0.4 


4.6±0.5 


6.4±1.5 


9±4 


0.25±0.06 


<0.82 


MM18423+5938 


34.1 ±8.8 


0.5 ±0.1 


0.79 ± 0.08 


11 ±3 


8±4 


0.41 ±0.10 


0.7 ±0.2 


SMMJ213511-0102 


22.4 ±0.6 


0.91 ±0.03 


0.95 ±0.01 


15.9 ±0.5 


14 ±6 


0.60 ±0.02 


0.37 ±0.02 


PSSJ2322+1944 


32.8±7.4 


0.71±0.11 


2.4±0.3 


4.9±1.0 


4.1±1.8 


0.19±0.04 


0.65±0.17 



Note. — All masses given in this table are corrected for lensing using fj, given in Table 4, Column 3. 
^ See Section 4.4. A temperature of 30 K has been assumed for all objects for which Tex could not be derived. 



Mcl = 5.706 X 10 



CI(l-O) 



[Mq] (4) 



where QiT^^) = 1 + 3e-^i/^- + Se-^^/T^^ jg the Ci 
partition function. ri = 23.6K and T2 = 62.5K are the 
energies above the ground stateQ- The above equation 
assumes optically thin C i emission and that both carbon 
lines are in LTE. 

Where available we use the excitation temperature de- 
rived in Section 4.4 above. For those sources where only 

12 In the case of J 1148+5251, for which only the Cl(2-1) line 
was observed, we use the corresponding equation given in WeiB et 
al. 2005b 



one carbon line has been measured we assume rox=30K. 
We note that the total neutral carbon mass is not a 
strong function of the assumed T^x unless the excitation 
temperature is below 20 K (unlikely to be the case in our 
targets, see also Figure 2 in Weifi et al. 2005a) so the 
exact choice of Tex is not critical for the derivation of C i 
masses. 

For completeness, we also derive molecular gas (H2) 
masses assuming the standard conversion factor for star- 
burst environments that should be applicable for SMGs 
and QSOs (aco=Mgas/i^o=0-8Mo(Kkms-i pc^)-!, 
Downes & Solomon 1998). We here assume that 
the J=3 line is slightly subthermally excited (i.e., 
^co(3-2)=0-9^^co(i-o)' ^^^^ 2005a). The 
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_^ CI{l-0) ^ C0(3-2) 



o 
a, 



,5^ 



10^ 10'" 10' 

log L'co(3-2) [K km pc2] 



O 



L pi/o_i\ L 



CI(2-1) ^C0(3-2) 



'^4 



10« 101" 
log L'co(3-2) [K km pc2] 



Fig. 13. — ^ci(i-o) ('^'^t) ^'^'^ ^ci(2-i) ('"'S'^*) ^ ^ function of ^co(3-2) ^'^'^ galaxies (stars indicate objects classified as quasars, 
circles: SMGs). The diagonal lines indicate ^ci~-'^cof3-2) ■'^Ci~^--'-^-'^cof3-2') 



CO luminosities (not corrected for magnification) are 
listed in Table 5, and the corresponding masses (here 
corrected for the magnification factors given in Table 4) 

in Tabic 6. Wc note that the molecular gas masses of 
the submillimeter galaxies would be higher by ~50% if 
we assumed the average iQQ^3_2)=0.6xi^Q^^_Qj derived 

by Ivison et al. (2010a). 

4.6. Atomic Carbon Abundances 

We can now derive the atomic carbon abundance rel- 
ative to Ha viaX[CI]/X[H2] = M(CI)/(6M(H2)). The 
carbon abundance derived this way is independent of the 
magnification (ignoring potentially differential magnifi- 
cation) and the applied cosmology. We find a carbon 
abundance of X[CI]/X[H2] « 8.4 ± 3.5 x lO--"^ for all 
sources (in agreement with earlier measurements by Ao 
et al., 2008, Weifi et al., 2005b). The carbon abundance 
for the SMGs would decrease if the 50% higher molecu- 
lar gas masses implied by Ivison et al. (2010a) were used. 
For comparison, Frerking et al. (1989) have derived an 
abundance of up to 2.2x10^"'' for high extinction regions 
in the Galactic star forming region Ophiuchus. Israel 
& Baas (2001, 2003) have estimated gas-phase abun- 
dances in nearby galaxies of order 4xl0~^ (their Table 
6). We conclude that the cold molecular gas (traced by 
C I and CO) is already substantially enriched in our sam- 
ple galaxies at high redshift. 

4.7. Cooling Contribution of Atomic Carbon 

As noted by Pety etal. (2004), a comparison of the 
C I and FIR luminosities shows that the cooling by Ci is 
negligible compared to the cooling by the dust continuum 
(average Lci/-^fir ~ 10~^). In Figure 14 (left) we plot 
the -^ci(i-o)/-^CO(3-2) ratio as a function of the cooling 
contribution of Ci (Lci/^'Fir)- We here chose to plot 



luminosity ratios to eliminate any uncertainties in mag- 
nification corrections (ignoring the possibility of differen- 
tial magnification). Sources referred to as SMGs appear 
to have higher ici(i_o)/^FiR I'atios than QSOs. How- 
ever, part of this offset is due to our choice of dust tem- 
peratures (Section 4.2.2): the lower Tjust in the SMGs 
(compared to the QSO sample) leads to lower LpiR, for a 
given S850;jm A^ix density, which in turn results in higher 
Lci(i_o)/^FiR ratios. To stick more to observables, we 
also plot the ratio of Lqi/ Lg^o^^ in the right panel of 
Figure 14 (£850pm=47rD^zyS850Mm; i^=352GHz). This is 
justified because for a given £fiRj -^850/jm is, to first or- 
der, independent of redshift for z>2 (e.g. Blain et al. 
2002). We still find an offset (higher Lci/Ls^Ofim) for the 
SMGs. This can be partly explained by a significant con- 
tribution of the central AGN to Lfir in the QSOs (e.g., 
discussion in Weifi et al., 2007a, for APM 08279-K5255, 
the source with the lowest ratio in our study). An al- 
ternative explanation is that the cooling contribution of 
Ci in the SMGs may be intrinsically higher than in the 
QSOs. 

4.8. Ci Ratios in the Context of PDR Models 

The PDR models of Gerin & Philhps (2000) and our 
average ratio Lci/LpiR ~ 10~^ suggest an average in- 
terstellar radiation field of order IO^xGq (where Go is 
the local Galactic interstellar ultraviolet radiation field) . 
Such high radiation fields have also been derived for other 
high-z sources (e.g., Maiolino et al. 2005, Ivison et al. 
2010b). 

Our derived average Ci(l-0)/Ci(2-l) line ratio (in 
units of L'} is 0.55±0.15 which translates to a ratio of 
2.7±0.9 in units of line intensity (which scales as i'^). For 
a 10*^x00 radiation field, the PDR models by Kaufman 
et al. (1999, their figure 8) and Meijerink et al. (2007, 
their figure 3) imply low H2 volume densities of order 
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Fig. 14. — Le/t; ^ci(i_()) /^co{3-2) ^ function of the cooling contribution of Cl (Lci/^fir)- Right: Same ^ci(i-o)/^CO(3-2) 

ratio, this time plotted as a function of ici/^850Mm to take out the FIR luminosity dependence on the assumed dust temperatures. SMGs 
are plotted as circles, QSOs as stars in both panels. 
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Fig. 15.— -f'ci(2_i)/-f'ci(i_o) 
nosity (corrected for lensing). 



ratio as a function of FIR lumi- 



10^ cm""^. This density is of similar order as the crit- 
ical density of both C i transitions which suggests that 
our derived excitation temperatures are, to first order, 
similar to the kinetic temperatures. Under this assump- 
tion, the average temperature of the C i lies in between 
the temperatures typically found in quiescent galactic re- 
gions ('Cirrus-like'), 15-20 K, and those found in active 
starforming environments (50-60 K). This could imply 
that the interstellar medium traced by our C i observa- 



tions is a mixture of a 'cold' and a 'warm' molecular gas 
component. We note however that the excitation tem- 
perture derived in equation 3 could be lower than the 
kinetic temperature if the C i excitation is not close to 
LTE. 

The C I line ratios are not in agreement with the pre- 
dictions by XDR models of Meijerink et al. (2007) which 
suggests that gas that is dominantly heated by an AGN 
should have a higher Ci(2-l)/Ci(l-0) ratio than ob- 
served (.^ci(2-i)/-^ci(i-o) >0-8, Fig. 3 in Meijerink). As 
shown in Figure 15 the C i ratio appears to be inde- 
pendent of the FIR luminosity and does not depend on 
the source type. Even in the case of APM 08279-1-5255, 
where the gas excitation is thought to be dominated by 
an AGN (Weifi et al. 2007), the C i(2-l)/Ci(l-0) ratio 
is lower than predicted for an XDR environment. 

5. DISCUSSION AND CONCLUSIONS 

We have presented new Ci observations of 10 high- 

redshift galaxies. All systems that have been looked at in 
C I emission to date (including studies in the literature) 
are very bright emitters both in the submillimeter regime 
and in CO line emission, and many of them are lensed. 
Even so, individual detections of the C i lines have proven 
to be challenging with current instrumentation. In total, 
we have targeted 16 C i lines in this study, and have 
detected 10 of them. Including these new data, there are 
now a total of 17 sources in the literature in which Ci 
has been observed (4 of which are not detected in neither 
the Ci(l-O) nor Ci(2 1) emission). 

Our main finding is that the C i properties of high- 
redshift systems do not differ significantly from what is 
found in low redshift systems, including the Milky Way. 
In addition, there are no major differences in C I proper- 
ties between the QSO and SMG-selected samples. We 



find that the Li 



CI(l-O)/ 



/L'qo ratios (0.29±0.12) are similar 



Atomic Carbon at High Redshift 



15 



to low-z galaxies (e.g., 0.2 ±0.2, Gerin & Phillips 2000). 

Measurements of both C i lines offer the unique op- 
portunity to constrain gas excitation tempertures of the 
molecular gas, independent of radiate transfer modeling. 
We derive an average carbon excitation temperature of 
29.1±6.3K for our sample. This temperature is lower 
than what is typically found in starforming regions in the 
local universe despite the fact that the sample galaxies 
have star formation rate surface densities on kpc scales 
of lOO's of M0 yr~^ kpc~^ (comparable to the most ex- 
treme starbursts in the local universe, but on larger spa- 
tial scales). However, the temperatures are roughly con- 
sistent with published dust temperatures of high-redshift 
starforming galaxies (Beelen et al. 2006, Kovacs et al., 
2006, 2010). Low carbon excitation as well as dust tem- 
peratures could indicate that the measurements include a 
significant amount of gas/dust unaffected by star forma- 
tion. However one should also keep in mind that the con- 
version of excitation temperature to kinetic temperature 
may be complicated by non-LTE excitation of atomic 
carbon. 

The C I abundances in our sample galaxies 
(X[Ci]/X[H2] = (8.4±3.5)xl0^5) are comparable, within 



the uncertainties, to what is found in local starforming 
environments (Section 4.6). This implies that the high-z 
galaxies studied here are significantly enriched in carbon 
on galactic scales, even though the look-back times are 
considerable (our average redshift of corresponds to 
an age of the universe of ^2 Gyr). 

In terms of the cooling budget the C i lines are a neg- 
ligible coolant (average Lci/-^^fir = (7.7 ±4.6) x 10~^). 
We find tentative evidence that this ratio may be ele- 
vated in the SMGs by a factor of a few compared to the 
QSOs, but a larger sample will be needed to beat down 
the current low-number statistics. The increase in sensi- 
tivity afforded by ALMA will be critical to increase the 
sample size and push C i studies at high redshift beyond 
the brightest systems that are accessible today. 



We thank the referee for comments that helped to im- 
prove the paper and thank Rowin Meijerink for useful 
discussions. This work is based on observations with the 
IRAM Plateau dc Bure Interferometer. IRAM is sup- 
ported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 



REFERENCES 



Alexander, D. M., Bauer, F. E., Chapman, S. C, Smail, I., Blain, 

A. W., Brandt, W. N., & Ivison, R. J. 2005a, ApJ, 632, 736 
Alexander, D. M., Smail, I., Bauer, F. E., Chapman, S. C, Blain, 

A. W., Brandt, W. N., & Ivison, R. J. 2005b, Nature, 434, 738 
Ao, Y., Weifi, A., Downes, D., Walter, F., Henkel, C, & Menten, 

K. M. 2008, A&A, 491, 747 
Barvainis, R., Maloney, P., Antonucci, R., & Alloin, D. 1997, 

ApJ, 484, 695 
Barvainis, R., & Ivison, R. 2002, ApJ, 571, 712 
Bayet, E., Gerin, M., Phillips, T. C, & Contursi, A. 2004, A&A, 

427, 45 

Beelen, A., Cox, P., Bonford, D. J., Dowell, C. D., Kovacs, A., 

Bertoldi, F., Omont, A., & Carilli, C. L. 2006, ApJ, 642, 694 
Benford, D. J., Cox, P., Omont, A., Phillips, T. C, & McMahon, 

R. G. 1999, ApJ, 518, L65 
Bertoldi, F., Cox, P., Neri, R., et al., 2003, A&A, 409, L47 
Blain, A. W., Smail, I., Ivison, R. J., Kneib, J. -P., & Frayer, 

D. T. 2002, Phys. Rep., 369, 111 
Blain, A. W., Barnard, V. E., & Chapman, S. C. 2003, MNRAS, 

338, 733 

Carilli, C. L., et al. 2002b, ApJ, 575, 145 
Carilli, C. L., et al. 2002a, AJ, 123, 1838 

Carilli, C. L., Lewis, G. F., Djorgovski, S. G., Mahabal, A., Cox, 

P., Bertoldi, F., & Omont, A. 2003, Science, 300, 773 
Casey, C. M., et al. 2009, MNRAS, 400, 670 
Chapman, S. C, et al. 2003, ApJ, 585, 57 
Cox, P., et al. 2002, A&A, 387, 406 

Daddi, E., Dannerbauer, H., Krips, M., Walter, F., Dickinson, M., 
Elbaz, D., & Morrison, G. E. 2009a, ApJ, 695, L176 

Daddi, E., et al. 2009b, ApJ, 694, 1517 

Danielson, A. L. R., et al. 2010, MNRAS, 1565 

Downes, D., Radford, J. E., Greve, A., Thum, C, Solomon, 
P. M., & Wink, J. E. 1992, ApJ, 398, L25 

Downes, D., Solomon, P. M., & Radford, S. J. E. 1995, ApJ, 453, 
L65 

Downes, D., & Solomon, P. M. 1998, ApJ, 507, 615 

Downes, D., Neri, R., Wiklind, T., Wilner, D. J., & Shaver, P. A. 

1999, ApJ, 513, LI 
Downes, D., & Solomon, P. M. 2003, ApJ, 582, 37 
Egami, E., Neugebauer, G., Soifer, B. T., Matthews, K., Ressler, 

M., Becklin, E. E., Murphy, T. W., Jr., & Dale, D. A. 2000, 

ApJ, 535, 561 

Fixsen, D.J., Bennett, C.L. & Mather, J.C. 1999, ApJ, 526, 207 
Frayer, D. T., Ivison, R. J., Scoville, N. Z., Yun, M., Evans, A. S., 

Smail, I., Blain, A. W., & Kneib, J.-P. 1998, ApJ, 506, L7 
Frayer, D. T., et al. 1999, ApJ, 514, L13 



Frerking, M. A., Keene, J., Blake, G. A., & Phillips, T. G. 1989, 

ApJ, 344, 311 
Gerin, M. & Phillips, T.G., 1998, ApJ 509, L14 
Gerin, M. & Phillips, T.G., 2000, ApJ 537, 644 
Genzel, R., Baker, A., Tacconi, L., Lutz, D., et al., 2003, ApJ, 

584, 633 

Greve, T. R., et al. 2005, MNRAS, 359, 1165 
Guilloteau, S., Omont, A., McMahon, R. G., Cox, P., & 

Petitjean, P. 1997, A&A, 328, LI 
Guilloteau, S., Omont, A., Cox, P., McMahon, R. G., & 

Petitjean, P. 1999, A&A, 349, 363 
Hainline, L. J., Scoville, N. Z., Yun, M. S., Hawkins, D. W., 

Frayer, D. T., & Isaak, K. G. 2004, ApJ, 609, 61 
Helou, G., Soifer, B. T., & Rowan-Robinson, M. 1985, ApJ, 298, 

L7 

Ikeda, M., Oka, T., Tatematsu, K., Sekimoto, Y., & Yamamoto, 

S. 2002, ApJS, 139, 467 
Isaak, K. G., Priddey, R. S., McMahon, R. G., Omont, A., 

Peroux, C, Sharp, R. G., & Withington, S. 2002, MNRAS, 

329, 149 

Israel, F. P., & Baas, F. 2001, A&A, 371, 433 
Israel, F.P, & Baas, F. 2002, A&A, 383, 82 
Israel, F. P., & Baas, F. 2003, A&A, 404, 495 
Israel, F.P 2005, Ap&SS, 295, 171 

Israel, F. P., Tilanus, R. P. J., & Baas, F. 2006, A&A, 445, 907 
Ivison, R. J., Smail, I., Le Borgne, J.-F., Blain, A. W., Kneib, 

J.-P., Bezecourt, J., Kerr, T. H., & Davies, J. K. 1998, 

MNRAS, 298, 583 
Ivison, R. J., Dunlop, J. S., Smail, I., Dey, A., Liu, M. C, & 

Graham, J. R. 2000, ApJ, 542, 27 
Ivison, R. J., Smail, I., Frayer, D. T., Kneib, J.-P., & Blain, 

A. W. 2001, ApJ, 561, L45 
Ivison, R. J., et al. 2002, MNRAS, 337, 1 

Ivison, R. J., Smail, I., Papadopoulos, P. P., Wold, I., Richard, J., 
Swinbank, A. M., Kneib, J.-P., & Owen, F. N. 2010b, MNRAS, 
404, 198 

Ivison, R. J., Papadopoulos, P. P., Smail, I., Greve, T. R., 
Thomson, A. P., Xilouris,E. M., & Chapman, S. C. 2010a, 
arXiv:1009.0749 

Ivison, R. J., et al. 2010c, A&A, 518, L35 

Kaufman, M. J., Wolfire, M. G., HoUenbach, D. J., & Luhman, 

M. L. 1999, ApJ, 527, 795 
Kneib, J.-P., Alloin, D., Mellier, Y., Guilloteau, S., Barvainis, R., 

& Antonucci, R. 1998, A&A, 329, 827 



16 



Walter, Weifi et al. 



Kneib, J. -P., van der Werf, P. P., Kraiberg Knudsen, K., Smail, I., 
Blain, A., Frayer, D., Barnard, V., & Ivison, R. 2004, MNRAS, 
349, 1211 

Kneib, J. -P., Neri, R., Smail, I., Blain, A., Sheth, K., van der 
Werf, P., & Knudsen, K. K. 2005, A&A, 434, 819 

Kovacs, A., Chapman, S. C, Dowell, C. D., Blain, A. W., Ivison, 
R. J., Smail, I., & Phillips, T. G. 2006, ApJ, 650, 592 

Kovacs, A., et al. 2010, ApJ, 717, 29 

Lestrade, J.-F., Combes, F., Salome, P., Omont, A., Bertoldi, F., 
Andre, P., & Schneider, N. 2010, A&A, 522, L4 

Lewis, G. F., Chapman, S. C, Ibata, R. A., Irwin, M. J., & 
Totten, E. J. 1998, ApJ, 505, LI 

Maiohno, R., et al. 2005, A&A, 440, L51 

Meijerink, R., Spaans, M., & Israel, F. P. 2007, A&A, 461, 793 
Menendez-Delmestre, K., et al. 2007, ApJ, 655, L65 
Menendez-Delmestre, K., et al. 2009, ApJ, 699, 667 
Neri, R., Genzel, R., Ivison, R.J et al. 2003, ApJ, 597, 113 
Ojha, R., Stark, A.A., Hsieh, H.H, et al. 2001, ApJ, 548, 253 
Oka, T., et al. 2005, ApJ, 623, 889 

Omont, A., Cox, P., Bertoldi, F., McMahon, R. G., CariUi, C, & 

Isaak, K. G. 2001, A&A, 374, 371 
Papadopoulos, P. P., Thi, W.-F., & Viti, S. 2004, MNRAS, 351, 

147 

Papadopoulos, P. P. 2005, ApJ, 623, 763 

Pety, J., Beelen, A., Cox, P., Downes, D., Omont, A., Bertoldi, 

F. , & CariUi, C. L. 2004, A&A, 428, L21 
Pope, A., et al. 2006, MNRAS, 370, 1185 

Riechers, D. A., Walter, F., Brewer, B. J., CariUi, C. L., Lewis, 

G. F., Bertoldi, F., & Cox, P. 2008b, ApJ, 686, 851 
Riechers, D. A., Walter, F., CariUi, C. L., Bertoldi , F., & 

Momjian, E. 2008a, ApJ, 686, L9 
Riechers, D. A., Walter, F., CariUi, C. L., & Lewis, G. F. 2009a, 

ApJ, 690, 463 
Riechers, D. A., et al. 2009b, ApJ, 703, 1338 

Robson, I., Priddey, R. S., Isaak, K. G., & McMahon, R. G. 2004, 

MNRAS, 351, L29 
Sanders, D. B., Mazzarella, J. M., Kim, D.-C, Surace, J. A., & 

Soifer, B. T. 2003, AJ, 126, 1607 
Schneider, N., Simon, R., Kramer, C, Kraemer, J., Stutzki, J., & 

Mookerjea, B. 2003, A&A, 406, 91 
Schoier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & 

Black, J. H. 2005, A&A, 432, 369 



Sheth, K., Blain, A. W., Kneib, J.-P., Frayer, D. T., van der 
Werf, P. P., & Knudsen, K. K. 2004, ApJ, 614, L5 

Solomon, P. M., Downes, D., & Radford, S. J. E. 1992, ApJ, 398, 
L29 

Solomon, P., Vanden Bout, P., 2005, ARA&A, 43, 677 

Spergel, D. N., et al. 2003, ApJS, 148, 175 

Spergel, D. N., et al. 2007, ApJS, 170, 377 

Stutzki, J., et al. 1997, ApJ, 477, L33 

Swinbank, A. M., et al. 2005, MNRAS, 359, 401 

Swinbank, A. M., et al. 2010, Nature, 464, 733 

Tacconi, L.J., Neri, R., Chapman, S.C. et al. 2006, ApJ, 640, 228 

Tacconi, L. J., et al. 2008, ApJ, 680, 246 

Takata, T., Sekiguchi, K., Smail, I., Chapman, S. C, Geach, J. E., 
Swinbank, A. M., Blain, A., & Ivison, R. J. 2006, ApJ, 651, 713 
Tielens, A. G. G. M., & Hollenbach, D. 1985a, ApJ, 291, 747 
Tielens, A. G. G. M., & Hollenbach, D. 1985b, ApJ, 291, 722 
Van der Werf, P. P., et al. 2010, arXiv: 1005.2877 
Venturini, S., & Solomon, P. M. 2003, ApJ, 590, 740 
Vernet, J., & Cimatti, A. 2001, A&A, 380, 409 
Villar-Martm, M., Fosbury, R. A. E., Binette, L., Tadhunter, 

C. N., & Rocca-Volmerange, B. 1999, A&A, 351, 47 
Wagg, J., Wilner, D. J., Neri, R., Downes, D., & Wiklind, T. 

2006, ApJ, 651, 46 
Walter, F., Weifi, A., & Scoville, N. 2002, ApJ, 580, L21 
Walter, F., Bertoldi, F., CariUi, C.L., et al. 2003, Nature, 424, 406 
Walter, F., CariUi, C, Bertoldi, F., Menten, K., Cox, P., Lo, 

K. Y., Fan, X., & Strauss, M. A. 2004, ApJ, 615, L17 
Weii3, A., Henkel, C, Downes, D., & Walter, F. 2003, A&A, 409, 

L41 

Weifi A., Downes, D., Henkel, C, Walter, F., 2005b, A&A, 429, 
L25 

Weifi A., Downes, D., Walter, F., Henkel, C. 2005a, A&A, 440, 
L45 

Weifi A., Downes, D., Neri, R., et al. 2007a, A&A 467, 955 
Weifi, A., Downes, D., Walter, F., & Henkel, C. 2007b, From 

Z-Machines to ALMA: (Sub)Millimeter Spectroscopy of 

Galaxies, 375, 25 
Weifi A., et al., 2011, in prep. 

White, G.J., Ellison, B., Claude, S., Dent, W.R.F. & Matheson, 

D. N., 1994, A&A 284, L23 



